An atrial fibrillation-associated regulatory region modulates cardiac Tbx5 levels and arrhythmia susceptibility.

Heart development and rhythm control are highly Tbx5 dosage-sensitive. TBX5 haploinsufficiency causes congenital conduction disorders, whereas increased expression levels of TBX5 in human heart samples has been associated with atrial fibrillation (AF). We deleted the conserved mouse orthologues of two independent AF-associated genomic regions in the Tbx5 locus, one intronic (RE(int)) and one downstream (RE(down)) of Tbx5. In both lines, we observed a modest (30%) increase of Tbx5 in the postnatal atria. To gain insight into the effects of slight dosage increase in vivo, we investigated the atrial transcriptional, epigenetic and electrophysiological properties of both lines. Increased atrial Tbx5 expression was associated with induction of genes involved in development, ion transport and conduction, with increased susceptibility to atrial arrhythmias, and increased action potential duration of atrial cardiomyocytes. We identified an AF-associated variant in the human RE(int) that increases its transcriptional activity. Expression of the AF-associated transcription factor Prrx1 was induced in Tbx5RE(int)KO cardiomyocytes. We found that some of the transcriptional and functional changes in the atria caused by increased Tbx5 expression were normalized when reducing cardiac Prrx1 expression in Tbx5RE(int)KO mice, indicating an interaction between these two AF genes. We conclude that modest increases in expression of dose-dependent transcription factors, caused by common regulatory variants, significantly impact on the cardiac gene regulatory network and disease susceptibility.

Patient-Specific TBX5-G125R Variant Induces Profound Transcriptional Deregulation and Atrial Dysfunction.

BACKGROUND: The pathogenic missense variant p.G125R in TBX5 (T-box transcription factor 5) causes Holt-Oram syndrome (also known as hand-heart syndrome) and early onset of atrial fibrillation. Revealing how an altered key developmental transcription factor modulates cardiac physiology in vivo will provide unique insights into the mechanisms underlying atrial fibrillation in these patients. METHODS: We analyzed ECGs of an extended family pedigree of Holt-Oram syndrome patients. Next, we introduced the TBX5-p.G125R variant in the mouse genome (Tbx5G125R) and performed electrophysiologic analyses (ECG, optical mapping, patch clamp, intracellular calcium measurements), transcriptomics (single-nuclei and tissue RNA sequencing), and epigenetic profiling (assay for transposase-accessible chromatin using sequencing, H3K27ac [histone H3 lysine 27 acetylation] CUT&RUN [cleavage under targets and release under nuclease sequencing]). RESULTS: We discovered high incidence of atrial extra systoles and atrioventricular conduction disturbances in Holt-Oram syndrome patients. Tbx5G125R/+ mice were morphologically unaffected and displayed variable RR intervals, atrial extra systoles, and susceptibility to atrial fibrillation, reminiscent of TBX5-p.G125R patients. Atrial conduction velocity was not affected but systolic and diastolic intracellular calcium concentrations were decreased and action potentials were prolonged in isolated cardiomyocytes of Tbx5G125R/+ mice compared with controls. Transcriptional profiling of atria revealed the most profound transcriptional changes in cardiomyocytes versus other cell types, and identified over a thousand coding and noncoding transcripts that were differentially expressed. Epigenetic profiling uncovered thousands of TBX5-p.G125R-sensitive, putative regulatory elements (including enhancers) that gained accessibility in atrial cardiomyocytes. The majority of sites with increased accessibility were occupied by Tbx5. The small group of sites with reduced accessibility was enriched for DNA-binding motifs of members of the SP (specificity protein) and KLF (Krüppel-like factor) families of transcription factors. These data show that Tbx5-p.G125R induces changes in regulatory element activity, alters transcriptional regulation, and changes cardiomyocyte behavior, possibly caused by altered DNA binding and cooperativity properties. CONCLUSIONS: Our data reveal that a disease-causing missense variant in TBX5 induces profound changes in the atrial transcriptional regulatory network and epigenetic state in vivo, leading to arrhythmia reminiscent of those seen in human TBX5-p.G125R variant carriers.

3D-cardiomics: A spatial transcriptional atlas of the mammalian heart.

Understanding the spatial gene expression and regulation in the heart is key to uncovering its developmental and physiological processes, during homeostasis and disease. Numerous techniques exist to gain gene expression and regulation information in organs such as the heart, but few utilize intuitive true-to-life three-dimensional representations to analyze and visualise results. Here we combined transcriptomics with 3D-modelling to interrogate spatial gene expression in the mammalian heart. For this, we microdissected and sequenced transcriptome-wide 18 anatomical sections of the adult mouse heart. Our study has unveiled known and novel genes that display complex spatial expression in the heart sub-compartments. We have also created 3D-cardiomics, an interface for spatial transcriptome analysis and visualization that allows the easy exploration of these data in a 3D model of the heart. 3D-cardiomics is accessible from http://3d-cardiomics.erc.monash.edu/.

Common Genetic Variants Contribute to Risk of Transposition of the Great Arteries.

RATIONALE: Dextro-transposition of the great arteries (D-TGA) is a severe congenital heart defect which affects approximately 1 in 4,000 live births. While there are several reports of D-TGA patients with rare variants in individual genes, the majority of D-TGA cases remain genetically elusive. Familial recurrence patterns and the observation that most cases with D-TGA are sporadic suggest a polygenic inheritance for the disorder, yet this remains unexplored. OBJECTIVE: We sought to study the role of common single nucleotide polymorphisms (SNPs) in risk for D-TGA. METHODS AND RESULTS: We conducted a genome-wide association study in an international set of 1,237 patients with D-TGA and identified a genome-wide significant susceptibility locus on chromosome 3p14.3, which was subsequently replicated in an independent case-control set (rs56219800, meta-analysis P=8.6x10-10, OR=0.69 per C allele). SNP-based heritability analysis showed that 25% of variance in susceptibility to D-TGA may be explained by common variants. A genome-wide polygenic risk score derived from the discovery set was significantly associated to D-TGA in the replication set (P=4x10-5). The genome-wide significant locus (3p14.3) co-localizes with a putative regulatory element that interacts with the promoter of WNT5A, which encodes the Wnt Family Member 5A protein known for its role in cardiac development in mice. We show that this element drives reporter gene activity in the developing heart of mice and zebrafish and is bound by the developmental transcription factor TBX20. We further demonstrate that TBX20 attenuates Wnt5a expression levels in the developing mouse heart. CONCLUSIONS: This work provides support for a polygenic architecture in D-TGA and identifies a susceptibility locus on chromosome 3p14.3 near WNT5A. Genomic and functional data support a causal role of WNT5A at the locus.

Neuregulin-1 enhances cell-cycle activity, delays cardiac fibrosis, and improves cardiac performance in rat pups with right ventricular pressure load.

OBJECTIVES: Right ventricular (RV) failure is a leading cause of death in patients with congenital heart disease. RV failure is kept at bay during childhood. Limited proliferation of cardiomyocytes is present in the postnatal heart. We propose that cardiomyocyte proliferation improves RV adaptation to pressure load (PL). We studied adaptation in response to increased RV PL and the role of increased cardiomyocyte cell cycle activity (CCA) in rat pups growing into adulthood. METHODS: We induced RV PL at day of weaning in rats (3 weeks; 30-40 g) by pulmonary artery banding and followed rats into adulthood (300 g). We performed histological analyses and RNA sequencing analysis. To study the effects of increased cardiomyocyte cell cycle activity, we administered neuregulin-1 (NRG1), a growth factor involved in cardiac development. RESULTS: PL induced an increase in CCA, with subsequent decline of CCA (sham/PL at 4 weeks: 0.14%/0.83%; P = .04 and 8 weeks: 0.00%/0.00%; P = .484) and cardiac function (cardiac index: control/PL 4 weeks: 4.41/3.29; P = .468 and 8 weeks: 3.57/1.44; P = .024). RNA sequencing analysis revealed delayed maturation and increased CCA pathways. NRG1 stimulated CCA (PL vehicle/NRG1 at 2 weeks: 0.62%/2.28%; P = .003), improved cardiac function (cardiac index control vs vehicle/NRG1 at 2 weeks: 4.21 vs 3.07/4.17; P = .009/.705) and postponed fibrosis (control vs vehicle/NRG1 at 4 weeks: 1.66 vs 4.82%/2.97%; P = .009/.078) in RV PL rats during childhood. CONCLUSIONS: RV PL during growth induces a transient CCA increase. Further CCA stimulation improves cardiac function and delays fibrosis. This proof-of-concept study shows that stimulation of CCA can improve RV adaptation to PL in the postnatal developing heart and might provide a new approach to preserve RV function in patients with congenital heart disease.

Epigenetic State Changes Underlie Metabolic Switch in Mouse Post-Infarction Border Zone Cardiomyocytes.

Myocardial infarction causes ventricular muscle loss and formation of scar tissue. The surviving myocardium in the border zone, located adjacent to the infarct, undergoes profound changes in function, structure and composition. How and to what extent these changes of border zone cardiomyocytes are regulated epigenetically is not fully understood. Here, we obtained transcriptomes of PCM-1-sorted mouse cardiomyocyte nuclei of healthy left ventricle and 7 days post myocardial infarction border zone tissue. We validated previously observed downregulation of genes involved in fatty acid metabolism, oxidative phosphorylation and mitochondrial function in border zone-derived cardiomyocytes, and observed a modest induction of genes involved in glycolysis, including Slc2a1 (Glut1) and Pfkp. To gain insight into the underlying epigenetic regulatory mechanisms, we performed H3K27ac profiling of healthy and border zone cardiomyocyte nuclei. We confirmed the switch from Mef2- to AP-1 chromatin association in border zone cardiomyocytes, and observed, in addition, an enrichment of PPAR/RXR binding motifs in the sites with reduced H3K27ac signal. We detected downregulation and accompanying epigenetic state changes at several key PPAR target genes including Ppargc1a (PGC-1α), Cpt2, Ech1, Fabpc3 and Vldrl in border zone cardiomyocytes. These data indicate that changes in epigenetic state and gene regulation underlie the maintained metabolic switch in border zone cardiomyocytes.

A Variant Noncoding Region Regulates Prrx1 and Predisposes to Atrial Arrhythmias.

[Figure: see text].

Variant Intronic Enhancer Controls SCN10A-short Expression and Heart Conduction.

BACKGROUND: Genetic variants in SCN10A, encoding the neuronal voltage-gated sodium channel NaV1.8, are strongly associated with atrial fibrillation, Brugada syndrome, cardiac conduction velocities, and heart rate. The cardiac function of SCN10A has not been resolved, however, and diverging mechanisms have been proposed. Here, we investigated the cardiac expression of SCN10A and the function of a variant-sensitive intronic enhancer previously linked to the regulation of SCN5A, encoding the major essential cardiac sodium channel NaV1.5. METHODS: The expression of SCN10A was investigated in mouse and human hearts. With the use of CRISPR/Cas9 genome editing, the mouse intronic enhancer was disrupted, and mutant mice were characterized by transcriptomic and electrophysiological analyses. The association of genetic variants at SCN5A-SCN10A enhancer regions and gene expression were evaluated by genome-wide association studies single-nucleotide polymorphism mapping and expression quantitative trait loci analysis. RESULTS: We found that cardiomyocytes of the atria, sinoatrial node, and ventricular conduction system express a short transcript comprising the last 7 exons of the gene (Scn10a-short). Transcription occurs from an intronic enhancer-promoter complex, whereas full-length Scn10a transcript was undetectable in the human and mouse heart. Expression quantitative trait loci analysis revealed that the genetic variants in linkage disequilibrium with genetic variant rs6801957 in the intronic enhancer associate with SCN10A transcript levels in the heart. Genetic modification of the enhancer in the mouse genome led to reduced cardiac Scn10a-short expression in atria and ventricles, reduced cardiac sodium current in atrial cardiomyocytes, atrial conduction slowing and arrhythmia, whereas the expression of Scn5a, the presumed enhancer target gene, remained unaffected. In patch-clamp transfection experiments, expression of Scn10a-short-encoded NaV1.8-short increased NaV1.5-mediated sodium current. We propose that noncoding genetic variation modulates transcriptional regulation of Scn10a-short in cardiomyocytes that impacts NaV1.5-mediated sodium current and heart rhythm. CONCLUSIONS: Genetic variants in and around SCN10A modulate enhancer function and expression of a cardiac-specific SCN10A-short transcript. We propose that noncoding genetic variation modulates transcriptional regulation of a functional C-terminal portion of NaV1.8 in cardiomyocytes that impacts on NaV1.5 function, cardiac conduction velocities, and arrhythmia susceptibility.

Early Postnatal Cardiac Stress Does Not Influence Ventricular Cardiomyocyte Cell-Cycle Withdrawal.

Congenital heart disease (CHD) is the most common birth defect. After birth, patients with CHD may suffer from cardiac stress resulting from abnormal loading conditions. However, it is not known how this cardiac burden influences postnatal development and adaptation of the ventricles. To study the transcriptional and cell-cycle response of neonatal cardiomyocytes to cardiac stress, we used a genetic mouse model that develops left ventricular volume overload within 2 weeks after birth. The increased volume load caused upregulation of the cardiac stress marker Nppa in the left ventricle and interventricular septum as early as 12 days after birth. Transcriptome analysis revealed that cardiac stress induced the expression of cell-cycle genes. This did not influence postnatal cell-cycle withdrawal of cardiomyocytes and other cell types in the ventricles as measured by Ki-67 immunostaining.

Genetic Dissection of a Super Enhancer Controlling the Nppa-Nppb Cluster in the Heart.

RATIONALE: ANP (atrial natriuretic peptide) and BNP (B-type natriuretic peptide), encoded by the clustered genes Nppa and Nppb, are important prognostic, diagnostic, and therapeutic proteins in cardiac disease. The spatiotemporal expression pattern and stress-induction of the Nppa and Nppb are tightly regulated, possibly involving their coregulation by an evolutionary conserved enhancer cluster. OBJECTIVE: To explore the physiological functions of the enhancer cluster and elucidate the genomic mechanism underlying Nppa-Nppb coregulation in vivo. METHODS AND RESULTS: By analyzing epigenetic data we uncovered an enhancer cluster with super enhancer characteristics upstream of Nppb. Using CRISPR/Cas9 genome editing, the enhancer cluster or parts thereof, Nppb and flanking regions or the entire genomic block spanning Nppa-Nppb, respectively, were deleted from the mouse genome. The impact on gene regulation and phenotype of the respective mouse lines was investigated by transcriptomic, epigenomic, and phenotypic analyses. The enhancer cluster was essential for prenatal and postnatal ventricular expression of Nppa and Nppb but not of any other gene. Enhancer cluster-deficient mice showed enlarged hearts before and after birth, similar to Nppa-Nppb compound knockout mice we generated. Analysis of the other deletion alleles indicated the enhancer cluster engages the promoters of Nppa and Nppb in a competitive rather than a cooperative mode, resulting in increased Nppa expression when Nppb and flanking sequences were deleted. The enhancer cluster maintained its active epigenetic state and selectivity when its target genes are absent. In enhancer cluster-deficient animals, Nppa was induced but remained low in the postmyocardial infarction border zone and in the hypertrophic ventricle, involving regulatory sequences proximal to Nppa. CONCLUSIONS: Coordinated ventricular expression of Nppa and Nppb is controlled in a competitive manner by a shared super enhancer, which is also required to augment stress-induced expression and to prevent premature hypertrophy.

Genome-Wide Analysis Identifies an Essential Human TBX3 Pacemaker Enhancer.

RATIONALE: The development and function of the pacemaker cardiomyocytes of the sinoatrial node (SAN), the leading pacemaker of the heart, are tightly controlled by a conserved network of transcription factors, including TBX3 (T-box transcription factor 3), ISL1 (ISL LIM homeobox 1), and SHOX2 (short stature homeobox 2). Yet, the regulatory DNA elements (REs) controlling target gene expression in the SAN pacemaker cells have remained undefined. OBJECTIVE: Identification of the regulatory landscape of human SAN-like pacemaker cells and functional assessment of SAN-specific REs potentially involved in pacemaker cell gene regulation. METHODS AND RESULTS: We performed Assay for Transposase-Accessible Chromatin using sequencing on human pluripotent stem cell-derived SAN-like pacemaker cells and ventricle-like cells and identified thousands of putative REs specific for either human cell type. We validated pacemaker cell-specific elements in the SHOX2 and TBX3 loci. CRISPR-mediated homozygous deletion of the mouse ortholog of a noncoding region with candidate pacemaker-specific REs in the SHOX2 locus resulted in selective loss of Shox2 expression from the developing SAN and embryonic lethality. Putative pacemaker-specific REs were identified up to 1 Mbp upstream of TBX3 in a region close to MED13L harboring variants associated with heart rate recovery after exercise. The orthologous region was deleted in mice, which resulted in selective loss of expression of Tbx3 from the SAN and (cardiac) ganglia and in neonatal lethality. Expression of Tbx3 was maintained in other tissues including the atrioventricular conduction system, lungs, and liver. Heterozygous adult mice showed increased SAN recovery times after pacing. The human REs harboring the associated variants robustly drove expression in the SAN of transgenic mouse embryos. CONCLUSIONS: We provided a genome-wide collection of candidate human pacemaker-specific REs, including the loci of SHOX2, TBX3, and ISL1, and identified a link between human genetic variants influencing heart rate recovery after exercise and a variant RE with highly conserved function, driving SAN expression of TBX3.

T-box transcription factor 3 governs a transcriptional program for the function of the mouse atrioventricular conduction system.

Genome-wide association studies have identified noncoding variants near TBX3 that are associated with PR interval and QRS duration, suggesting that subtle changes in TBX3 expression affect atrioventricular conduction system function. To explore whether and to what extent the atrioventricular conduction system is affected by Tbx3 dose reduction, we first characterized electrophysiological properties and morphology of heterozygous Tbx3 mutant (Tbx3+/-) mouse hearts. We found PR interval shortening and prolonged QRS duration, as well as atrioventricular bundle hypoplasia after birth in heterozygous mice. The atrioventricular node size was unaffected. Transcriptomic analysis of atrioventricular nodes isolated by laser capture microdissection revealed hundreds of deregulated genes in Tbx3+/- mutants. Notably, Tbx3+/- atrioventricular nodes showed increased expression of working myocardial gene programs (mitochondrial and metabolic processes, muscle contractility) and reduced expression of pacemaker gene programs (neuronal, Wnt signaling, calcium/ion channel activity). By integrating chromatin accessibility profiles (ATAC sequencing) of atrioventricular tissue and other epigenetic data, we identified Tbx3-dependent atrioventricular regulatory DNA elements (REs) on a genome-wide scale. We used transgenic reporter assays to determine the functionality of candidate REs near Ryr2, an up-regulated chamber-enriched gene, and in Cacna1g, a down-regulated conduction system-specific gene. Using genome editing to delete candidate REs, we showed that a strong intronic bipartite RE selectively governs Cacna1g expression in the conduction system in vivo. Our data provide insights into the multifactorial Tbx3-dependent transcriptional network that regulates the structure and function of the cardiac conduction system, which may underlie the differences in PR duration and QRS interval between individuals carrying variants in the TBX3 locus.

Trait-associated noncoding variant regions affect TBX3 regulation and cardiac conduction.

Genome-wide association studies have implicated common genomic variants in the gene desert upstream of TBX3 in cardiac conduction velocity. Whether these noncoding variants affect expression of TBX3 or neighboring genes and how they affect cardiac conduction is not understood. Here, we use high-throughput STARR-seq to test the entire 1.3 Mb human and mouse TBX3 locus, including two cardiac conduction-associated variant regions, for regulatory function. We identified multiple accessible and functional regulatory DNA elements that harbor variants affecting their activity. Both variant regions drove gene expression in the cardiac conduction tissue in transgenic reporter mice. Genomic deletion from the mouse genome of one of the regions caused increased cardiac expression of only Tbx3, PR interval shortening and increased QRS duration. Combined, our findings address the mechanistic link between trait-associated variants in the gene desert, TBX3 regulation and cardiac conduction.

Identification of Functional Variant Enhancers Associated With Atrial Fibrillation.

RATIONALE: Genome-wide association studies have identified a large number of common variants (single-nucleotide polymorphisms) associated with atrial fibrillation (AF). These variants are located mainly in noncoding regions of the genome and likely include variants that modulate the function of transcriptional regulatory elements (REs) such as enhancers. However, the actual REs modulated by variants and the target genes of such REs remain to be identified. Thus, the biological mechanisms by which genetic variation promotes AF has thus far remained largely unexplored. OBJECTIVE: To identify REs in genome-wide association study loci that are influenced by AF-associated variants. METHODS AND RESULTS: We screened 2.45 Mbp of human genomic DNA containing 12 strongly AF-associated loci for RE activity using self-transcribing active regulatory region sequencing and a recently generated monoclonal line of conditionally immortalized rat atrial myocytes. We identified 444 potential REs, 55 of which contain AF-associated variants (P<10-8). Subsequently, using an adaptation of the self-transcribing active regulatory region sequencing approach, we identified 24 variant REs with allele-specific regulatory activity. By mining available chromatin conformation data, the possible target genes of these REs were mapped. To define the physiological function and target genes of such REs, we deleted the orthologue of an RE containing noncoding variants in the Hcn4 (potassium/sodium hyperpolarization-activated cyclic nucleotide-gated channel 4) locus of the mouse genome. Mice heterozygous for the RE deletion showed bradycardia, sinus node dysfunction, and selective loss of Hcn4 expression. CONCLUSIONS: We have identified REs at multiple genetic loci for AF and found that loss of an RE at the HCN4 locus results in sinus node dysfunction and reduced gene expression. Our approach can be broadly applied to facilitate the identification of human disease-relevant REs and target genes at cardiovascular genome-wide association studies loci.

Identification of atrial fibrillation associated genes and functional non-coding variants.

Disease-associated genetic variants that lie in non-coding regions found by genome-wide association studies are thought to alter the functionality of transcription regulatory elements and target gene expression. To uncover causal genetic variants, variant regulatory elements and their target genes, here we cross-reference human transcriptomic, epigenomic and chromatin conformation datasets. Of 104 genetic variant regions associated with atrial fibrillation candidate target genes are prioritized. We optimize EMERGE enhancer prediction and use accessible chromatin profiles of human atrial cardiomyocytes to more accurately predict cardiac regulatory elements and identify hundreds of sub-threshold variants that co-localize with regulatory elements. Removal of mouse homologues of atrial fibrillation-associated regions in vivo uncovers a distal regulatory region involved in Gja1 (Cx43) expression. Our analyses provide a shortlist of genes likely affected by atrial fibrillation-associated variants and provide variant regulatory elements in each region that link genetic variation and target gene regulation, helping to focus future investigations.

Conserved NPPB+ Border Zone Switches From MEF2- to AP-1-Driven Gene Program.

BACKGROUND: Surviving cells in the postinfarction border zone are subjected to intense fluctuations of their microenvironment. Recently, border zone cardiomyocytes have been specifically implicated in cardiac regeneration. Here, we defined their unique transcriptional and regulatory properties, and comprehensively validated new molecular markers, including Nppb, encoding B-type natriuretic peptide, after infarction. METHODS: Transgenic reporter mice were used to identify the Nppb-positive border zone after myocardial infarction. Transcriptome analysis of remote, border, and infarct zones and of purified cardiomyocyte nuclei was performed using RNA-sequencing. Top candidate genes displaying border zone spatial specificity were histologically validated in ischemic human hearts. Mice in which Nppb was deleted by genome editing were subjected to myocardial infarction. Chromatin accessibility landscapes of border zone and control cardiomyocyte nuclei were assessed by using assay for transposase-accessible chromatin using sequencing. RESULTS: We identified the border zone as a spatially confined region transcriptionally distinct from the remote myocardium. The transcriptional response of the border zone was much stronger than that of the remote ventricular wall, involving acute downregulation of mitochondrial oxidative phosphorylation, fatty acid metabolism, calcium handling, and sarcomere function, and the activation of a stress-response program. Analysis of infarcted human hearts revealed that the transcriptionally discrete border zone is conserved in humans, and led to the identification of novel conserved border zone markers including NPPB, ANKRD1, DES, UCHL1, JUN, and FOXP1. Homozygous Nppb mutant mice developed acute and lethal heart failure after myocardial infarction, indicating that B-type natriuretic peptide is required to preserve postinfarct heart function. Assay for transposase-accessible chromatin using sequencing revealed thousands of cardiomyocyte lineage-specific MEF2-occupied regulatory elements that lost accessibility in the border zone. Putative injury-responsive enhancers that gained accessibility were highly associated with AP-1 (activator protein 1) binding sites. Nuclear c-Jun, a component of AP-1, was observed specifically in border zone cardiomyocytes. CONCLUSIONS: Cardiomyocytes in a discrete zone bordering the infarct switch from a MEF2-driven homeostatic lineage-specific to an AP-1-driven injury-induced gene expression program. This program is conserved between mouse and human, and includes Nppb expression, which is required to prevent acute heart failure after infarction.

Transcriptome analysis of mouse and human sinoatrial node cells reveals a conserved genetic program.

The rate of contraction of the heart relies on proper development and function of the sinoatrial node, which consists of a small heterogeneous cell population, including Tbx3+ pacemaker cells. Here, we have isolated and characterized the Tbx3+ cells from Tbx3+/Venus knock-in mice. We studied electrophysiological parameters during development and found that Venus-labeled cells are genuine Tbx3+ pacemaker cells. We analyzed the transcriptomes of late fetal FACS-purified Tbx3+ sinoatrial nodal cells and Nppb-Katushka+ atrial and ventricular chamber cardiomyocytes, and identified a sinoatrial node-enriched gene program, including key nodal transcription factors, BMP signaling and Smoc2, the disruption of which in mice did not affect heart rhythm. We also obtained the transcriptomes of the sinoatrial node region, including pacemaker and other cell types, and right atrium of human fetuses, and found a gene program including TBX3, SHOX2, ISL1 and HOX family members, and BMP and NOTCH signaling components conserved between human and mouse. We conclude that a conserved gene program characterizes the sinoatrial node region and that the Tbx3+/Venus allele provides a reliable tool for visualizing the sinoatrial node, and studying its development and function.

Specialized impulse conduction pathway in the alligator heart.

Mammals and birds have a specialized cardiac atrioventricular conduction system enabling rapid activation of both ventricles. This system may have evolved together with high heart rates to support their endothermic state (warm-bloodedness) and is seemingly lacking in ectothermic vertebrates from which first mammals then birds independently evolved. Here, we studied the conduction system in crocodiles (Alligator mississippiensis), the only ectothermic vertebrates with a full ventricular septum. We identified homologues of mammalian conduction system markers (Tbx3-Tbx5, Scn5a, Gja5, Nppa-Nppb) and show the presence of a functional atrioventricular bundle. The ventricular Purkinje network, however, was absent and slow ventricular conduction relied on trabecular myocardium, as it does in other ectothermic vertebrates. We propose the evolution of the atrioventricular bundle followed full ventricular septum formation prior to the development of high heart rates and endothermy. In contrast, the evolution of the ventricular Purkinje network is strongly associated with high heart rates and endothermy.

EMERGE: a flexible modelling framework to predict genomic regulatory elements from genomic signatures.

Regulatory DNA elements, short genomic segments that regulate gene expression, have been implicated in developmental disorders and human disease. Despite this clinical urgency, only a small fraction of the regulatory DNA repertoire has been confirmed through reporter gene assays. The overall success rate of functional validation of candidate regulatory elements is low. Moreover, the number and diversity of datasets from which putative regulatory elements can be identified is large and rapidly increasing. We generated a flexible and user-friendly tool to integrate the information from different types of genomic datasets, e.g. ATAC-seq, ChIP-seq, conservation, aiming to increase the ease and success rate of functional prediction. To this end, we developed the EMERGE program that merges all datasets that the user considers informative and uses a logistic regression framework, based on validated functional elements, to set optimal weights to these datasets. ROC curve analysis shows that a combination of datasets leads to improved prediction of tissue-specific enhancers in human, mouse and Drosophila genomes. Functional assays based on this prediction can be expected to have substantially higher success rates. The resulting integrated signal for prediction of functional elements can be plotted in a build-in genome browser or exported for further analysis.

OccuPeak: ChIP-Seq peak calling based on internal background modelling.

UNLABELLED: ChIP-seq has become a major tool for the genome-wide identification of transcription factor binding or histone modification sites. Most peak-calling algorithms require input control datasets to model the occurrence of background reads to account for local sequencing and GC bias. However, the GC-content of reads in Input-seq datasets deviates significantly from that in ChIP-seq datasets. Moreover, we observed that a commonly used peak calling program performed equally well when the use of a simulated uniform background set was compared to an Input-seq dataset. This contradicts the assumption that input control datasets are necessary to fatefully reflect the background read distribution. Because the GC-content of the abundant single reads in ChIP-seq datasets is similar to those of randomly sampled regions we designed a peak-calling algorithm with a background model based on overlapping single reads. The application, OccuPeak, uses the abundant low frequency tags present in each ChIP-seq dataset to model the background, thereby avoiding the need for additional datasets. Analysis of the performance of OccuPeak showed robust model parameters. Its measure of peak significance, the excess ratio, is only dependent on the tag density of a peak and the global noise levels. Compared to the commonly used peak-calling applications MACS and CisGenome, OccuPeak had the highest sensitivity in an enhancer identification benchmark test, and performed similar in an overlap tests of transcription factor occupation with DNase I hypersensitive sites and H3K27ac sites. Moreover, peaks called by OccuPeak were significantly enriched with cardiac disease-associated SNPs. OccuPeak runs as a standalone application and does not require extensive tweaking of parameters, making its use straightforward and user friendly. AVAILABILITY: http://occupeak.hfrc.nl.